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The effects of long-range anisotropie elastic deformations on electronic structure in supercon- 
ductors are analyzed within the framework of the Bogohubov-de Gennes equations. Cases of twin 
boundaries and isolated defects are considered as illustrations. We find that the superconducting 
order parameter is depressed in the regions where pronounced lattice deformation occurs. The calcu- 
lated locai density of states suggests that the electronic structure is strongly modulated in response 
to lattice deformations, and propagates to longer distances. In particular, this allows the trapping 
of low-lying quasiparticle states around defects. Some of our predictions can be directly tested by 
STM experiments. 
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In many complex electronic materials such as cuprates, 
manganites, ferroelastic martensites, and titanates, un- 
expected and puzzling multiscale modulations of charge, 
spin, polarization, and strain variables bave been re- 
vealed by high resolution microscopy 0. The nonuni- 
form textures found in these doped materials indicate 
that their origin is intrinsic: they arise from coupling 
between various degrees of freedom. The textures fun- 
damentally affect locai and mesoscopic electronic, niag- 
netic and structural properties, which are centrai to the 
functionality of correlated electronic materials. There is 
ampie evidence for significant coupling amongst the elec- 
tronic degrees of freedom with the lattice distortion in 
cuprates and manganites. The charge carrier doping can 
act as a locai stress to deform surrounding unit cells 
We might employ a Landau- Ginzburg (LG) theory to 
study the coupling between the electronic (Cooper pair) 
and lattice (strain tensor) degrees of freedom in super- 
conductors. However, the LG theory can only describe 
the long wavelength behavior. The locai electronic prop- 
erties and lattice distortion necessitate a treatment at 
the atomic scale. Recently, we have developed an atomic 
scale theory for determining lattice distortions by using 
strain related variables and their constraint equations • 
This now enables a systematic study of the influence of 
strain on electronic wavefunctions. Hcrc we apply a mi- 
croscopie theory to study the order parameter and locai 
quasiparticle properties in both s and d-wave supercon- 
ductors. 

We consider the foUowing model on a square lattice: 
^ = - 5Z *y + ^(e, - lj)c\^c,a 
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Here Cia annihilates an electron of spin cr on site i. The 
quantities Ci and fi are the on-site impurity potential 
(if any) and the chcmical potential, respectivcly. The 
hopping integrai Uj is modified by the lattice distor- 
tion. The electron-lattice coupling is approximated by 



t, 



whcrc t?. 



is the bare hopping integrai. 



Cij is the lattice-distortion variable, and a is the coupling 
Constant. In our nearest neighbor realization, the bare 
hopping integrai t^^ is t for ncarcst neighbor sites and zero 
otherwisc. Specifically, we take the form of the lattice dis- 
tortion to be: €ij = [|(Rj-fdj) — (Ri-|-di)|/|Rj— Ri| — 1], 
where {Ri} are the undistorted lattice coordinates and 
{di} the lattice displacement vectors with respect to 
{Ri}. We assume an effective superconducting gap func- 



tion given by Ai 



2 



{ci^Cji—CiiCj-f), where Uij = U5i. 



(i.e., attractive Hubbard-C/ model) for s-wave supcrcon- 
ductivity and Uij = V5i-\.^j (with 7 specifying the near- 
est ncighbors to the i-th site) for d-wavc superconduc- 
tivity. By performing a Bogoliubov-Valatin transforma- 
tion, we may diagonalize our Hamiltonian by solving the 
Bogoliubov-dc Gennes (BdG) cquation 0: 



E 



(2) 



subject to the self-consistency conditions for the super- 
conducting (SC) order parameter (OP): 



tanh 



/ En 
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Here the single particle Hamiltonian reads Uij = —tij + 
{ci — n)Sij. The quasiparticle wavefunction, correspond- 
ing to the eigenvalue i5„, consists of electron (w") and 
hole (w") amplitudes. The quasiparticle energy is mea- 
sured with respect to the chemical potential. 

We solve the BdG equations self-consistcntly by start- 
ing off with an initial gap function. After exactly diago- 
nalizing Eq. the obtained wavefunction is substituted 
info Eq. Q to computo a ncw gap function. We then 
use this as an input to repeat the abovc procedure until 
a desired convergence is achieved. Below, we report our 
results for two types of locai lattice distortions at zero 
temperature — a superlattice formed by twin boundaries 
and a single defect. We measure the length and energy 
in units of ao (the undistorted lattice Constant) and t. 
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FIG. 1: Strain-eg mode for a periodic twinned microstruc- 
ture (a) and a single defect (c) together with their correspond- 
ing displacement configurations [(b) and (d)] within the high- 
lighted window. Nl=2,2x 32. 



The chemical potential /i = and no extrinsic impurity 
scattering is introduced (e^ = 0). The pairing interaction 
for both the s-wave {U) and d-wave {V) superconductors 
is taken to be 3. The typical system size is Nl ~ 32 x 32 
and periodic boundary conditions are applied. When 
the locai quasiparticle density of states (LDOS) is com- 
puted, we implement a much larger system using the 
above small system as a supercell. In the absence of dis- 
tortions, Eq. Q is readily solved by resorting to trans- 
lational invariance: the resulting quasiparticle energy 
E-k = Vsk+^k 'with = — 2(cosfca. + cosky) — and 
the energy gap Ak = A^q or Ak = ■^^^(cosfc^ — cosfc.y) 
for the s or d-wa.ve superconductor respectively. For 
the given parameter values, we obtain Ago = 0.85 and 
A^o = l-'''. In both cases, the superconducting coher- 
ence length is ~ 2ao. These values are exaggerated when 
compared to real materials: The choice of values is mo- 
tivated by the desire to enhance the effect of the lattice 
distortions. 

Before proceeding to the BdG calculation sketched 
above, we generate (following two atomic scale lat- 
tice distortions. These distortions arise from long-range 
anisotropie interactions between strains to maintain the 
compatibility constraints. The es (square to rectangle) 
strain mode for a periodic twinned microscopie struc- 
turc and the corrcsponding atomic displaccmcnts are de- 
pictcd in Fig. ^a-b). The domain consists of rcctangu- 
lar distortion (rcd and blue) separated by a domain wall 
whcre 63 = 0. Note that similar but more "realistic" mi- 
crostructures bave also been obtaincd using Monte Carlo 
techniques 0. Similar quantities for a single defect are 
displayed in pancls (c-d). Near the defect, four alternat- 
ing distortions are formed in (dover leaf) "c?-wave like" 
pattern. 

In Fig. 13 we show the spatial variation of the SC OP 
induced by the deformation of Fig. ^a) in both s and 
d-wave superconductors. In both types of superconduc- 
tors, the OP is lowered within the domain and is ele- 
vated at the domain wall (Fig. |3a-b)). The magnitude 
of the OP is depressed in comparison to an undistorted 
square lattice since the lattice deformation changes the 
band structure, leading to a reduction in normal den- 
sity of states at the Fermi energy. Even at the domain 
wall, where the strain induced deformation is weakest, 



FIG. 2: Spatial variation of the SC OP for periodic twin 
boundaries displayed in Fig. H^a) — (a) The s-wave OP in an 
s-wave superconductor, and (b) the d-wave and (c) extended 
s-wave components of the OP in a d-wave superconductor. 
The electron-lattice coupling Constant a = 3. 



the amplitude of the enhanced OP is smaller than its 
value in an undistorted square lattice. This is due to the 
confincmcnt from the two neighboring domains. In an 
s-wave superconductor. the relative orbitai motion be- 
tween two paired electrons has an angular momentum 
l — 0, which has the highcst symmctry and gencratcs no 
subdominant OPs with lower symmetries. However, for a 
d-wave superconductor, the relative motion between two 
paired electrons has the angular momentum 1 = 2, which 
upon scattering from any inhomogeneity can generate a 
subdominant OP with symmetry not lower than that of 
l = 2. The d-wavc OP is defined by Ad{i) = (Aì^ì+ì, -I- 
Ai.i_£ — Ai_i^y — Ai^i_j;)/4, whereas the extended s-wave 
OP is As{i) = (Aj,,+£ + Aj,,_i + Ai^,+y + Aj,,_{;)/4. In 
a twinned domain of a d-wave SC, a subdominant ex- 
tended s-wave componcnt is gencratcd in a real combi- 
nation d ± s. Because the symmetries of two twinned 
domains are reflected info each other with respect to the 
twin boundary, the relative phase between the d- and s- 
wave components switches by tt when a twin boundary 
is crossed (Fig. Eie)). It has been argued phenomeno- 
logically 01 that a locai time-reversal-symmetry-breaking 
state exists at a twin or grain boundary of YBa2Cu307_5. 
Within numcrical accuracy, our rcsult shows a real ad- 
mixturc of the d-wavc and s-wave components of the OP. 
A d+is pairing state was also found at the {110}-orientcd 
surface or interface of a d-wave superconductor in early 
work 0. A cruciai difference between the twin bound- 
ary of Fig. [3 and the {110}-oriented surface of earlier 
work is that the dominant d-wave component reaches a 
maximum at twin boundaries, whereas it is strongiy sup- 
pressed at the {110}-oriented surface or interface. Ex- 
perimentally, the existence of a time-reversal-symmetry- 
brcaking pairing state in high- Te cupratcs is the subject 
of current debate ■ 

As another example, we show in Fig.|3|the spatial vari- 
ation of the superconducting OP around the single defect 
(Fig.^c)) in both the s-wave and d-wave superconductor 
cases. The OP is depressed at the center of the defect, 
and reaches its defect-free bulk value at the length scale 
^0- Notice that for a lattice-deformation defect, which 
affects the locai electron hopping integrai, the OP has 
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FIG. 3: Spatial variation of the SC OP for a single defect dis- 
played in Fig. 0^c) — (a) The s-wave OP in an s-wave super- 
conductor, and (b) the d-wave and (c) extended s-wave com- 
ponents of the OP in a d-wave superconductor. The electron- 
lattice coupling Constant a = 3. 



a minimum at four sitcs surrounding the dcfcct center. 
It is difFerent from the case of an externaUy substituted 
unitary impurity. where the minimum of OP is located 
only at the impurity site itself 01. The range of influ- 
ence of such a defect can be very large depending on 
the strength of electron-lattice coupling — the elasticity 
propagates the electronic response. The d-wave energy 
gap has a sign change at the nodal directions of the es- 
sentially cylindrical Fermi surface, but the d-wave OP 
does not exhibit such a sign change in real space. Whcn 
the defect is introduced, an extended s-wave component 
of the OP is induced whcn the dominant d-wave compo- 
nent is dcprcssed at the defect. Strikingly, the induced 
s-wave component has a sign change across the diagonals 
of the square lattice, i.e., sgn[cos(2é')], where is the az- 
imuthal angle with respect to the crystalline x axis. This 
is a direct manifestation of the d-wave pairing symme- 
try in real space. This fcaturc can be understood from 
a phenomenological LG free-energy density functional 
with a two-component SC OP: J- = as|Asp -|- adjA^P + 
/3i|A,|4 + /32|A,|4-H/33|A.|2|Arf|2-)-/34(Af A^-f-Af A2) + 
7« I VA, |2 -h7d| VA<j|2 4-7.d[9^ A,a^ A* - A,9j, A* -hc.c] , 
where we take a, to be always positive while ad = 
ctdo{T/Tdo — 1) such that there exists only a single tran- 
sition into a d-wave pairing state in a homogeneous sys- 
tem. Whcn a dcfcct is introduced, the d-wave component 
is depressed. Since the s-wave component itself should 
be very small, the term A*^ A^ + A^^A^ is only a higher- 
order correction. Therefore, it is the mixed-gradient term 
that induces the s-wave component and also determines 
the relative phase to be or tt. The sign change is also 
evident by exchanging the x and y components of the 
position coordinate in the mixed-gradient term. 

Once the self-consistcncy for the order parameter is 
obtaincd. we calculatc the LDOS: 



p,(i?) = -^[K|V'(i?-i?«) + k"lV'(i? + £^n)] , (4) 



where f'{E) is the derivative of the Fermi distribution 
function with respect to the energy. The LDOS deter- 
mines the differential tunneling conductance, measurable 
by STM experiments [13 ■ Figure shows the LDOS at 
a domain wall for both types of superconductors, where 
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FIG. 4; The locai density of states at a twin boundary in s- 
wave (a) and d-wave superconductors. Also shown the LDOS 
(black lines) for a single domain. The electron-lattice coupling 
Constant a — Z. 



the modulation of the superconducting OP forms a su- 
perlattice, with its maximum at the domain wall play- 
ing the role of an off-diagonal potential barrier (A^j in 
Eq. l|2Jl)- For an s-wave superconductor, the quasiparti- 
cles are gapped away with their energy below the mini- 
mum SC OP. Outside the minimum of the pair potential, 
energy bands are formed by the quasiparticle scattering 
off the off-diagonal energy barriers at the domain walls. 
Except for the gap about the Fermi energy (i? = 0), this 
is reminiscent of the electronic structure in semiconduc- 
tor supcrlattices (e.g., alternating GaAs/GaAlAs layers). 
Interestingly, the bottom of the oscillation pattern fol- 
lows the LDOS (black line) of a system formed by a single 
rectangular domain. Similar oscillations are obtained for 
the d-wave superconductor. However, the bottom of the 
oscillations do not foUow the single domain DOS (black 
line). In addition, weak subgap peaks (labeled by arrows 
in Fig. 0^b)) appear symmetrically in the LDOS on the 
domain wall but are absent in the single-domain LDOS. 
We speculate that these resonant states are due to the 
gradient of the s-wave gap component induced inside the 
domain. 

As shown in Fig. we have also calculated the LDOS 
near the center of a single defect. The depression of the 
SC OP at the defect makes a quantum-well-like profile 
of the energy gap. The size and depth of the well is de- 
termined by the electron-lattice coupling Constant. Be- 
cause of the difference between the s-wave and d-wave 
pairings, the s-wave potential well is closed everywhere, 
whereas the d-wave well has four slits along the diagonals 
of the square lattice. In the s-wave superconductor, the 
well is shallow and small for weak coupling, which can- 
not trap low-lying quasiparticle bound states; for strong 
coupling constants, the wcU is deep and large so that 
subgap quasiparticle bound states are induced (the red 
and blue Hnes of Fig. E^a) and (c)). The energy of these 
low-lying states must be inbetween the bottom and edge 
of the well. Therefore, it is notable that the energy of 
these subgap states is shifted toward the Fermi surface 
as the electron-lattice coupling is increased (the blue line 
in Fig. |SIa) and (c)). The situation hcrc is also diffcrcnt 



3 



LA-UR-03-3177 



Aprii 29, 2003 



Submittcd to Physical Rcvicw Lettera 




E 



FIG. 5: The locai density of states near the center of a defect 
in s-wave (left column) and d-wave superconductors (right 
column). The distance of the measured point away from the 
defect is labeled by its coordinate. The electron-lattice cou- 
pling constants are a = 3 (red lines) and 10 (blue lines). Also 
shown is the defect-free LDOS (black lines). 



from an s-wave vortex core, where the OP at the core 
center must vanish such that the low-lying bound states 
are always trapped The electronic structure at the 
defect in a d-wave superconductor becomes even richer: 
For a = 3 (weak coupHng as compared to the band width 
of the non-deformed square lattice), the lattice distortion 
plays the role of a weak defect for the quasiparticle scat- 
tering. In this resonant peak with a dip exactly 
at the Fermi energy is seen (the red line in Fig. E{b)). 
The overall peak comes from the scattering of quasipar- 
ticles off the single-particle off-diagonal potential (i.e., 
locai change of the hopping integrai as a response to the 
lattice deformation). This lattice-dcformation induced 
resonance state also exhibits Friedel oscillations. Typi- 
cally, the peak structure appears in the LDOS at (0,0) 
(We label the four sites surrounding the defect center 
by (0,0), (1,0), (1,1), (0,1)) and (-2,-2). In contrast to 
the case of an extrinsic on-site potential-scattering im- 
purity [ll| . the LDOS spectrum is symmctric since even 
the locai particlc-holc symmetry is prcscrved here. For 
a = 10 (strong coupling), the d-wave OP is almost fuUy 
depressed (less than 0.03), while the maximum of the in- 
duced s-wave OP reachcs 0.085. The locai ofF-diagonal 
potential becomes more finite ranged, which causcs a 
locai change of the band width. The 'resonant' peaks 
are pushed to higher energies (~ ±0.3) (the blue line of 
Fig. Od)). Furthermore, small shoulders appear dose to 
the Fermi energy (the blue hnes of Fig. |Sfb) and (d)), 
which are precursors of new Andreev resonance states. 
We have also computed the LDOS without imposing self- 
consistcncy on the OP and found that the doublc-pcak 
structure is V-shapcd with no cxistcncc of the shoulders. 



This leads us to speculate that the new Andreev reso- 
nance states come from the confinement of the induced 
s-wave OP. However, these states are stili delocalized be- 
causG the quasiparticles can leak out of the well through 
the slits along the diagonal directions where the induced 
s-wave component vanishes. Ali these features are unique 
to an elastic defect in a d-wave superconductor with short 
coherence length. 

In conclusion, we studied the effects of elastic lattice 
deformation on the nanoscale electronic structure in su- 
perconductors. We have shown that the SC OP is de- 
pressed in the regions where the lattice deformation ex- 
ists. The calculated LDOS suggests that the electronic 
structure is strongly modulatcd in responso to the lattice 
deformation. In particular, it is possiblc to trap low-lying 
quasiparticle states around the dcfccts. Imagcs of these 
states will manifcst the undcrlying long-rangc anisotropie 
lattice deformation. These prcdictions can bc dircctly 
tested by STM experiments in new functional supercon- 
ducting materials. Our approach is readily extended to 
other elastic textures and SC symmetries. Self-consistent 
coupling of the elastic and SC textures on an equal foot- 
ing will be pursued elsewhere. It would also be inter- 
esting to study the electronic response by using Monte 
Carlo gencratcd microstructurcs as input, which consti- 
tutes future work. 
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